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Dynamic scaling analyses are performed in the spin-glass phase of the ± J Ising, the XY, and 
the Heisenberg models in three dimensions. We found a crossover from the critical dynamics 
to the ground-state dynamics in the Ising model and the Heisenberg model. The ground-state 
dynamics of the Ising model is characterized by an activation law with a finite energy gap: 
the typical time diverges exponentially. On the other hand, the typical time in the Heisenberg 
model diverges algebraically with the inverse temperature. Algebraic relaxation with a finite 
dynamic exponent is observed after the typical time in both models. The ground-state dynamic 
exponent is estimated to be zo — 13, which is common to both models. There is no crossover 
in the XY model. The critical dynamics is considered to continue to the ground-state. 
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Randomness and frustration produce nontrivial phe- 
nomena in spin glasses. 1 Theoretical investigations in 
spin glasses are mostly concentrated into two topics. One 
is an issue regarding a picture of the spin-glass (SG) 
phase. There have been many arguments whether it is ex- 
plained by the replica-symmetry-breaking(RSB) theory 2 
or it is explained by the droplet theory. 3 The second topic 
is an argument of existence or nonexistence of a spin- 
glass transition in short-range Edwards- Anderson(EA) 
models. There is a consensus that the SG transition oc- 
curs in the Ising model in three dimensions. 4,5 However, 
the issue is not settled on the vector spin models. A chi- 
rality mechanism insists that a chiral-glass(CG) transi- 
tion occurs without the SG transition. 6, 7 However, there 
are several investigations suggesting a simultaneous tran- 
sition. 8-11 

In our recent paper, 12 the CG and the SG transition 
temperatures, T sg and T cg , are resolved to be different 
in the Heisenberg model. We have performed a dynamic 
scaling analysis of the glass susceptibilities. Particularly, 
we paid attention to a temperature dependence of the 
dynamic exponent, z(T). Since the exponent is a key 
variable in the dynamic scaling hypothesis, a careful de- 
termination of z(T) increases accuracy of the scaling re- 
sults. 

It is known that the dynamic exponent is proportional 
to the inverse temperature in the spin-glass phase. 13-15 

I T 

The behavior is also observed experimentally. 16 It has 
been considered an outcome of the slow dynamics in the 
spin-glass phase. However, it is found recently 1 that the 
temperature dependence also appears in the paramag- 
netic phase and there is no anomaly at the transition 
temperature. It is certified to be a general feature of the 
spin-glass models in three dimensions. 12 

It is important to investigate whether or not the tem- 
perature dependence of the dynamic exponent continues 
to the zero temperature. If it continues, a dynamics of 



the ground state is smoothly connected to the param- 
agnetic state. We must consider a theory of the spin- 
glass dynamics including both the ground state and the 
paramagnetic state. If it does not continue and there 
is a crossover temperature, then we may consider the 
higher and the lower temperature regions separately. The 
ground-state (or close to the ground-state) dynamics is 
intrinsically different from the critical dynamics near the 
spin-glass transition temperature. Hukushima et al. 18 re- 
ported such a crossover in the four-dimensional ± J Ising 
model. 

In this paper, we study the low-temperature phase of 
the ± J Ising, the XY, and the Heisenberg models in three 
dimensions. A dynamic scaling analysis is performed on 
the spin-glass correlation length. The temperature de- 
pendences of the dynamic exponent are obtained. We 
find a dynamic crossover in the Ising model and in the 
Heisenberg model. The spin dynamics is essentially dif- 
ferent above and below the crossover temperature, T cr . 

We consider an EA model on a simple cubic lattice. 
The Hamiltonian is written as 



T~i — ^ JijSi 

(hj) 



(2) 



The interactions, Jij, take the two values 



J and — J 

with the same probability. The temperature, T, is scaled 
by J. Linear lattice size is denoted by L. A total num- 
ber of spins is N = L x L x (L + l), and skewed periodic 
boundary conditions are imposed. The spins are updated 
by the single-spin-flip algorithm using the Metropolis 
probability. 

We start simulations with random spin configurations 
(T = oo). The temperature is quenched to a finite value 
at the first Monte Carlo step. Physical quantities we ob- 
serve are the spin-glass (chiral-glass) correlation length, 
which is estimated from the correlation functions by a 
single exponential fitting. The relaxation functions are 
averaged over independent Monte Carlo simulations. 

A typical lattice size in the Ising simulations is 71. A 
replica number is 128, and typical numbers of random 
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Fig. 1. (Color online) (a) Nonequilibrium relaxation functions of 
the spin-glass correlation length in the Ising model. The temper- 
atures of data are denoted aside the fitting lines. A relaxation 
function at T = is plotted by circle symbols, (b) The dynamic 
scaling plot. When T > T cr = 0.5, data are plotted with symbols 
against (t/4) T / T =g . Below T cr , they are plotted with lines against 
10 4 t/exp[4/T]. 



bond samples are more than one hundred. In the XY 
simulations, the lattice size is 55, and the replica num- 
ber is 64. The sample numbers are mostly around one 
hundred. In the Heisenberg simulations, the lattice size 
is 59, and the replica number is 32 or 64. The sample 
numbers are between 10 and 100. 

Figure 1(a) shows relaxation functions of the spin-glass 
correlation length in the Ising model. They show alge- 
braically diverging behaviors. We estimate the nonequi- 
librium dynamic exponent, z, from the logarithmic slope, 
since 



m = £o x (t/r) 1 ^ 



(3) 



from the dynamic scaling relation. The dynamic expo- 
nent smoothly depends on the temperature when T > 
T cr = 0.5. 

The relaxation function changes the behavior at T cr . 
It shows a plateau below this temperature. As the tem- 
perature decreases, the plateau range becomes longer. A 
spin state is considered to be trapped in a local energy 
minimum. After some time steps, it eventually escapes 
from the valley and shows an algebraic divergence, whose 
dynamic exponent is independent from the temperature. 
This is a distinct difference from the higher temperature 
region. A temperature effect only appears in the dynam- 
ics of escaping from the valley. A dynamics after the es- 
cape is independent from the temperature. It is natural 
to consider that this temperature-independent dynam- 
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Fig. 2. (Color online) Dynamic scaling plots of the SG and the 
CG correlation lengths in the XY model. CG data are scaled 
only below T sg . 



ics continues to the ground state. Therefore, this is a 
crossover from the critical dynamics around the transi- 
tion temperature to the ground state dynamics. 

We perform a dynamic scaling analysis on the spin- 
glass correlation length. A basic scaling form is the 
dynamic scaling relation, Eq. (3). When T > T cr , a 
typical time, r, is a constant and the dynamic expo- 
nent, z, is proportional to the inverse temperature. The 
scaling analysis is performed by plotting £(t) against 
(t/r) T / T ^. 19 ' 20 We use an estimate of T sg = 1.18(1) 
obtained in our previous paper. 12 Here, r is an only 
scaling parameter. From the logarithmic slope we obtain 
the dynamic critical exponent at the transition temper- 
ature, z sg . When T < T cr , a typical time, r, depends 
on the temperature and the dynamic exponent, z, is a 
constant value. We suppose an activation law for r as 
r = Toexp[A/T], where A is an average energy barrier 
of the valley. In this scaling, A is a scaling parameter. 

Figure 1(b) shows results of the scaling analysis. The 
high-tempcrature scaling is achieved by choosing t = 4. 
The results are plotted by symbols. Each temperature 
data fall on a single line. From the slope we obtain 
z sg = 5.1(1). 12 The low-temperature scaling is achieved 
by choosing A = 4J. Algebraic increases after the 
plateau ride on a single line. The dynamic exponent of 
this process is estimated to be zq = 13(1). 

A value of A = 4J is a special number. It is an en- 
ergy difference of a local single-spin flip when two spins 
are energetically favorable and four spins are energeti- 
cally unfavorable among the six nearest neighbors. This 
energy valley is due to a local single-spin flip. It is not 
a global excitation valley. The temperature only controls 
this local activation. If we renormalize a Monte Carlo 
step by the activation time, we can neglect a dynamic 
process of trapping by the local energy minimum. Then, 
the temperature-independent ground-state relaxation of 
the spin-glass correlation length is observed. It is an aver- 
aged global increase of the frozen domain. The estimated 
zq is the ground-state dynamic exponent. 

Figures 2 shows the dynamic scaling result of the SG 
and the CG correlation lengths in the XY model. A scal- 
ing of z(T) cx 1/T is supposed. We use estimates of tran- 
sition temperatures obtained in our previous paper. 12 
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Fig. 3. (Color online) Nonequilibrium relaxation functions of the 
SG(a) and the CG(b) correlation lengths in the Heisenberg 
model, (c) The dynamic scaling plots. When T > T cr = 0.07, 
data are plotted with symbols against (f /3000) T//Ts s . Below T cr , 
they are plotted with lines against t/T 1,7 . Algebraic divergence 
with z sg and that with zq are also plotted for guide for eyes. 



The scaling plot of CG is possible only below T sg . The 
correlation times of SG and CG are same. The dynamic 
critical exponent at each transition temperature is also 
consistent with each other. Therefore, the dynamic be- 
havior in the low-temperature phase is explained by a 
critical dynamics of the spin-glass transition. It is con- 
sidered to continue to the ground state. There is no 
crossover temperature. Of course, there is a possibility 
of a crossover below the lowest temperature we have in- 
vestigated. 

Figure 3 shows the results in the Heisenberg model. 
Relaxation functions of the correlation lengths show the 
same behavior as in the Ising model. The logarithmic 
slope depends on the temperature when T > T cr ~ 0.07. 
We have also performed a simulation at T = 0, which is 
just a spin-quenching dynamics. The correlation length 
shows algebraic divergence with zo = 13, which is the 



ground-state dynamic exponent. This value coincides 
with the Ising value. When T < T cr , the correlation 
length exhibits algebraic divergence with the ground- 
state dynamic exponent in the long-time limit. There is 
no qualitative difference between the spin and the chiral- 
ity 

A clear difference between the Ising model and the 
Heisenberg model is a lack of plateau in the low- 
temperature relaxation. The correlation length increases 
algebraically before it shows the ground-state relaxation. 
It is considered as a consequence of continuous distribu- 
tion of the energy barrier in the local spin flips. Typ- 
ical length scale where the ground-state relaxation ap- 
pears depends on the temperature. There is finite ther- 
mal fluctuation even in the spin-glass phase. It affects the 
short-range spin correlations. The ground-state behavior 
is only observed when the spin-glass correlation length 
exceeds this thermally-fluctuating domain size, which be- 
comes large as the temperature increases. In the Ising 
model, the domain size is considered to be constant due 
to a discrete energy barrier. 

The scenario mentioned above is well-depicted in the 
dynamic scaling result in Fig. 3(c). We have plotted 
the spin-glass correlation length against the renormal- 
ized Monte Carlo steps. When T > T cr , the scaling anal- 
ysis is same as the Ising model. Typical time, r = 3000, is 
constant and the dynamic exponent depends on the tem- 
perature as z(T) = 5.1 x T sg /T. A value of T sg = 0.18 
is taken from our previous paper. 12 When T < T cr , the 
relaxation function is scaled with the typical time which 
algebraically diverges as T~ 1,7 . It is also considered to be 
a consequence of continuous distribution of the energy 
barrier A. We observe the ground-state dynamics after 
renormalizing the Monte Carlo steps by this typical time. 
The convergence to the ground-state relaxation becomes 
later as the temperature increases. It is also noted that 
the relaxation before the convergence is the critical relax- 
ation with the dynamic exponent z sg . Therefore, what is 
observed in this dynamic scaling is a crossover from the 
critical dynamics in the short-time regime to the ground- 
state dynamics in the long-time regime. 

Figure 4 summarizes the temperature dependence of 
the dynamic exponent. The Ising model and the Heisen- 
berg model exhibit the same behavior. The value of the 
ground-state dynamic exponent agrees with each other. 
The similarity suggests that the symmetry of both 
models is responsible to the dynamic crossover behavior. 
We may consider that there is no crossover in the XY 
model due to a lack of this symmetry. 

To conclude, we have investigated the dynamic behav- 
ior in the spin-glass phase. We observed a crossover from 
the critical dynamics to the ground-state dynamics in 
the Ising model and in the Heisenberg model. The dy- 
namic exponent z(T) changes the temperature depen- 
dence. The algebraic ground-state dynamics is observed 
after renormalizing Monte Carlo steps by a typical time. 
The ground-state dynamic exponent is estimated to be 
zo ~ 13, which is common to the two models. The typ- 
ical time diverges exponentially in the Ising model, and 
it diverges algebraically in the Heisenberg model. There 
is no crossover in the XY model. The dynamic property 
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Fig. 4. (Color online) Temperature dependences of the nonequi- 
librium dynamic exponent, z(T), and the static exponent, "f/u, 
in the Ising model. 



of the low-temperature phase is explained solely by the 
critical dynamics. 

A motivation of this study is to answer a question 
in regard to the nature of the spin-glass phase: RSB or 
droplet. The present results alone cannot give an exclu- 
sive answer. However, it is possible to interpret our re- 
sults by two theories. In a viewpoint of the droplet theory, 
a crossover from the critical dynamics to the ground- 
state dynamics is observed. Even at very low temper- 
atures, the spin-glass correlation length diverges alge- 
braically with time. It may be the T = criticality of 
the droplet theory. 3 However, the ground-state dynamic 
exponent is estimated to be finite and is independent of 
the temperature. This result is different from the origi- 
nal droplet theory. 3 In the XY model, the dynamic ex- 
ponent seems to diverge at T — 0. However, it is difficult 
to ensure the criticality. A diverging dynamic exponent 
suggests that the SG/CG correlation length does not di- 



verge or at most logarithmically diverges. In a viewpoint 
of the RSB theory, the present simulation time is too 
short to observe the RSB behavior. Only an effect of a 
local energy valley is observed within our algorithm and 
computer facility. Up to the present, our results are bet- 
ter explained by the droplet theory. 

A dynamic crossover observed in this paper suggests 
that the spin-glass ground-state is accessible even at a 
finite temperature if below the crossover temperature. 
Particularly, the finite energy barrier in the Ising model 
may be overcome by a finite quantum fluctuation of the 
transverse field. It can be a reason why the quantum an- 
nealing method works well in the spin-glass problems. 21 
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